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Abstract 



in 



We observed CO J = 3 — 2 emission from the "water fountain" sources, which exhibit high-velocity col- 
limated stellar jets traced by H2O maser emission, with the Atacama Submillimeter Telescope Experiment 
(ASTE) 10 m telescope. We detected the CO emission from two sources, IRAS 16342-3814 and 
dr IRAS 18286-0959. The IRAS 16342-3814 CO emission exhibits a spectrum that is well fit to a Gaussian 

q \ profile, rather than to a parabolic profile, with a velocity width (FWHM) of 158 ±6 km s _1 and an inten- 

5h ' sity peak at Vlsr= 50±2 km s _1 . The mass loss rate of the star is estimated to be ~ 2.9 x 1O _5 M yr _1 . 

c/3 \ Our morpho-kinematic models suggest that the CO emission is optically thin and associated with a bipolar 

outflow rather than with a (cold and relatively small) torus. The IRAS 18286—0959 CO emission has a 
velocity width (FWHM) of 3.0 ±0.2 km s _1 , smaller than typically seen in AGB envelopes. The narrow 
velocity width of the CO emission suggests that it originates from either an interstellar molecular cloud or 
' a slowly- rotating circumstellar envelope that harbors the H2O maser source. 

G\ ! Key words: stars: AGB and post-AGB — stars: individual (IRAS 16342-3814, IRAS 18286-0959) 

O 

1. Introduction detection of water fountain sources. This implies that 

these jets should have just been launched at the final evo- 
The process in which mass outflow shapes planetary lutionary stage of the dying stars. Note that the H2O 
nebulae (PNe) is still unresolved and needs to be eluci- maser emission is excited when the bipolar tips of the 
■ dated. One of the striking features of some PNe is the stellar jet strike into the ambient CSE that was produced 
bipolarity in their morphological and kinematical struc- from the previous spherically symmetric stellar mass loss. 
v~j ■ tures, some of which clearly exhibit high velocity bipo- Therefore, in order to more reliably estimate the stellar 
KnJ I lar jets (Vj et >100 km s _1 ). Such bipolar jets have mass loss rate, the dynamical age, and the whole morpho- 
been found at the beginning of the post asymptotic gi- logical and kinematical structure of the jet, a mapping 
ant branch (AGB) phase and even at the end of the AGB observation of thermal emission such as CO emission is 
phase (e.g. Sahai & Trauger 1998; Imai et al. 2002; crucial. 

Imai 2007). In particular, a small fraction of stellar jets The first detection of CO emission from a water 
have velocities higher than a typical expansion velocity of fountain source was reported with the Arizona Radio 
a circumstellar envelope (CSE) found in 1612 MHz OH Observatory 10 m telescope towards IRAS 16342-3814 
maser emission (Vj ct >30 km s" 1 ), and are traced by H 2 (He et al. 2008). Usually it is difficult to detect such CO 
maser emission. Such high velocity jets are called "water emission towards the water fountain sources because they 
fountains" (Imai 2007). As of 2009, 13 such objects have are located close to the Galactic plane with heavy contam- 
been identified (Walsh et al. 2009; Suarez et al. 2009). ination from the interstellar CO emission or they are too 
Some spatio-kinematical structures of the jets have been distant (D >2 kpc). Here we present results of the first 
elucidated by very long baseline interferometry (VLBI) systematic CO J = 3-2 emission observations towards the 
observations of the H 2 mascrs (Imai et al. 2002; Imai water fountain sources with the Atacama Submillimeter 
et al. 2005; Imai et al. 2007; Imai 2007; Boboltz & Telescope Experiment (ASTE) 10 m telescope. In addi- 
Marvel 2007; Claussen et al. 2009). Interestingly, all tion, we also show results of CO J = 1 - emission obser- 
of the estimated apparent dynamical ages of the jets are vations towards W 43A and IRAS 18286-0959 with the 
about 100 yr or shorter, which is consistent with the rare Nobeyama 45 m telescope. In section 2 and 3, we describe 
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in detail the observations and results, respectively. In sec- 
tion 4, we discuss the mass loss rate and the morphology 
and kinematics of IRAS 16342-3814, whose CO J = 3 - 2 
emission is detected. 

2. Observations and data reduction 

2.1. ASTE observations 

The ASTE observations of 12 CO J = 3 - 2 emis- 
sion at 345.79599 GHz were made on June 20 and 
21, 2008 during LST 13:00-22:00. The FWHM beam 
size of the ASTE telescope is 22" at the observed fre- 
quency For antenna temperature calibration, the CO 
emission from IRC+10216 was observed in a 10 min 
scan on both days. Comparing the antenna tempera- 
ture of IRC+10216 with that obtained with the Caltech 
Submillimcter Observatory (CSO) 10 m telescope, we ob- 
tained an antenna aperture efficiency rj mb =0.59 with rea- 
sonable agreement with those found in previous ASTE 
observations. An antenna pointing check was made every 
^2 hr. The antenna pointing was confirmed to be stable 
within ~5". The system temperature was between 200 
and 400 K (single side band). The received signals were 
down-converted in frequency and transferred into three 
base band channels (BBCs), each of which has a band 
width of 512 MHz, corresponding to a velocity width of 
445 km s _1 at 345 GHz. The center frequencies of the 
BBCs were split by 100 MHz to check the intrinsic CO 
emission, which should be detected in different spectral 
channels but at the same velocity ranges in the all BBCs. 
We used the MAC spectrometer to obtain a spectrum with 
1024 spectral channels, corresponding to a frequency and 
velocity spacings of 500 kHz and 0.45 km s _1 , respec- 
tively. Reduction of the spectral data was made using 
the NewStar package developed at the Nobcyama Radio 
Observatory. 

Table 1 gives parameters of the water fountain sources 
and the observations. The observations were made in 
the antenna-position switching mode; the number of 
points observed on and around the target source (>1 for 
the cross-scan mode) is given in Column 10. Because 
IRAS 18286-0959, W 43A, and IRAS 18460-0151 are 
located very close to the Galactic plane with intense back- 
ground molecular emission, five-points or nine-points cross 
scans were adopted, in which the observed points were 
separated by 10"-20"(see sect. 3). Column 11 in table 1 
gives root-mean-square noise levels of the spectra. Note 
that the raw spectrum after integration had an emission- 
free baseline with a linear gradient. This situation was 
nearly ideal for detecting any velocity component of CO 
emission. The linear gradient baseline was subtracted to 
obtain the final spectrum. 

2.2. Nobeyama ^5 m observations 

The Nobcyama 45 m observations of 12 CO J = 1 — 
emission at 115.271204 GHz were made towards W 43A 
on December 1-2, 2008 during LST 16:00-21:00 and to- 
wards IRAS 18286-0959 on April 29, 2008 with a beam 
size of about 15". A single beam SIS receiver and the 



5x5 beam focal plane array receiver system (BEARS, 
Sunada et al. 2000) were used to make five point and 25 
point scans around the target sources, respectively. The 
grid spacing of BEARS was 41", while that of the five 
point scans was set to 30". The system temperatures were 
between 400 and 600 K (single side bands). To observe 
W 43A and IRAS 18286-0959, both of the two receiving 
systems and BEARS were used, respectively. For signals 
received with the SIS receiver and BEARS, the acousto- 
optic spectrometers (AOSs) and the digital spectrometers 
(Sorai et al. 2000) were used, respectively, to obtain spec- 
tra that have 2048 and 1024 spectral channels, yielding 
velocity spacings of 0.33 km s _1 and 0.16 km s _1 (in a 
two-channel binding mode) each at 115 GHz, respectively. 
In order to find intrinsic emission from the target sources, 
the spectra obtained from four points around each target 
were averaged and the average spectrum was subtracted 
from that obtained at the target position. 

3. Results 

We found CO J = 3 - 2 emission with ASTE 
towards IRAS 15445-5449, IRAS 15544-5332, 
IRAS 16342-3814, IRAS 18286-0959, W 43A, 
IRAS 18460-0151, and IRAS 18596+0315. From 
IRAS 15445-5449, IRAS 15544-5332, W 43A, 
IRAS 18460-0151, and IRAS 18596+0315, however, 
the emission must originate from background sources 
because no intrinsic CO emission is detected. If it 
indeed exists, it should be detected only at the exact 
position and systemic velocity of the source. It has 
been observationally found that the center velocities 
of H2O and 1612 MHz OH maser spectra in the water 
fountain sources are coincident with each other to within 
1-2 km s _1 (Imai et al. 2002; Imai 2007). This implies 
that the systemic velocities of the driving sources of 
these masers as well as the CO envelope/flow should 
coincide within the above range. The systemic velocities 
are cited from previous H2O and SiO maser spectra 
(see table 1). At positions around W 43A, the same 
velocity components with equal antenna temperature in 
each component are found. Separately, we have observed 
12 CO, 13 CO J = 1 - and HCN emission towards 
W 43A with the Nobeyama 45 m telescope, but obtained 
negative detection (with upper limits of 0.2 K, 0.1 K 
in antenna temperatures for 12 CO and 13 CO J = 1 — 
lines, respectively). The J = 3 — 2 CO emission towards 
IRAS18286-0959 and IRAS16342-3814 is detected at the 
precise position and velocity of each source. Towards 
other water fountain sources, no CO emission was 
detected in the position-switching observations. 

Fig. 1 shows the CO J = 3 — 2 spectrum towards 
IRAS 16342-3814 (hereafter abbreviated as 116342). To 
improve the sensitivity, three MAC spectra covering differ- 
ent velocity ranges are synthesized although the intrinsic 
CO emission is already identified in each of the spectra. 
The synthesized spectrum is well fit to a Gaussian pro- 
file with a velocity width (FWHM) of 158 + 6 km s" 1 
and a peak intensity of I pC ak = 32.9 + 1.0 mK at Vlsr= 
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50 ± 2 km s _1 with respect to the local standard of rest 
(LSR). Remarkably, the velocity width is larger than that 
of 1612 MHz OH maser emission (~120 km s _1 , Sahai 
et al. 1999; Claussen et al. 2009) and comparable to 
the full width of the H 2 maser spectrum (e.g., Likkel 
et al. 1992; Claussen et al. 2009). The origin of the wide 
wing of the spectrum is discussed in sect. 4.2. 

H2O masers in IRAS 18286—0959 were discovered by 
Deguchi et al. (2007). Negative detections of HCO+, 
Ha53, and NH3 lines toward this object suggest that they 
are not associated with any dense molecular cloud or com- 
pact HII region. Although a plenty of velocity compo- 
nents in the H2O masers may imply their association of 
a massive young stellar object, their spatio-kinematical 
structure clearly exhibits collimated jets as seen in water 
fountain sources (Imai 2007). Fig. 2 shows the result of 
the CO J = 3 - 2 observation towards IRAS 18286-0959. 
Towards IRAS 18286—0959 and nearby positions, five ve- 
locity components are found in the CO J = 3 — 2 spec- 
tra. Only the velocity component at Vlsr=65 km s _1 
has its antenna temperature maximum at the position of 
IRAS 18286-0959 (fig. 2a and b). An off-source spectrum 
was synthesized from the spectra at four positions around 
this source, and was subtracted from the on-source spec- 
trum. In the five velocity components, only one velocity 
component at Vlsr— 65 km s _1 was identified as possible 
intrinsic emission from IRAS 18286-0959 (Fig. 2c). As 
shown in figure 3, the same velocity component is also 
found in CO J = 1 — emission. However, this component 
has a velocity width (FWHM) of only 3 km s _1 (fig. 2d), 
which is smaller than those typically found in AGB en- 
velopes (~10 km s^ 1 , e.g. Kemper et al. 2003) and much 
smaller than the velocity range of the H2O maser emis- 
sion in the same source (Deguchi et al. 2007; Imai 2007). 
Our unpublished European VLBI Network (EVN) data 
shows that 1612 MHz OH maser emission is detected at 
Vlsr=39.5 km s" 1 towards IRAS 18286-0959. We be- 
lieve that the OH maser component is one of double- 
peaked components as typically seen in stellar 1612 MHz 
OH maser sources. If the CO and OH lines are generated 
from the common circumstellar envelope, a terminal ex- 
pansion velocity is derived to be ~25 km s _1 , which is 
inconsistent with the narrow line width of the CO emis- 
sion. Using Equation 1 in sect. 4.1, we derive a mass loss 
rate to be ~ 3.5 x 10~ 7 M Q yr -1 , inconsistent with the 
water fountain sources being at the final AGB or early 
post-AGB phase. We think that the CO emission may be 
associated with a background young stellar object seen in 
the same direction. Alternatively, it may be associated 
with a rotating circumstellar disk surrounding the cen- 
tral object that collimates the H 2 maser flow. In fact, 
the existence of such circumstellar disks has already been 
proposed (e.g., Jura & Kahanc 1999; Nakashima 2006). 
In order to unambiguously elucidate the origin of the CO 
J = 3 — 2 and J = 1 — lines, mapping observations with 
radio interferometers are indispensable. 



4. Discussion 

4.1. Stellar mass loss rate of IRAS 16342-3814 

In the present paper, we attempt to estimate a mass 
loss rate of the 116342 CO J = 3 — 2 outflow. In spite of 
the limited quality of the 13 CO J = 2 — 1 spectrum, He 
et al. (2008) suggests that the CO J = 2 — 1 line is optically 
thin on the basis of a low 12 CO antenna temperature. In 
addition, the CO J = 3 — 2 line is axisymmctric with re- 
spect to the systemic velocity. Thus it is expected that 
the CO J = 3 — 2 line is optically thin. In this case, using 
Eq. (6) and (8) of Knapp & Morris (1985) and a correc- 
tion for the beam size and the different CO transition (cf. 
Olofsson et al. 1993; Groenewegen et al. 1999), a mass 
loss rate is estimated in units of solar masses per year as 
follows, 



M gas = 4.55x10" 



Log(W/0.04)s(J) 
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fc^J 6 (DB) 
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Here Tmb is the antenna temperature of the CO emis- 
sion in Kelvin, V^p the expansion velocity of the CO 
emission in km s _1 , D the source distance in parsec, 
B the beam size of the telescope, /co the abundance 
of CO molecules relative to H 2 , s(J) a correction factor 
for J — ► J — 1 transition, and W the ratio of the 4.6 /xm 
flux to that emitted by a blackbody of temperature of 
2000 K and radius of 5xl0 13 cm (i.e., W=l). For 116342, 
Tmb =33 mK, D =2 kpc, and B =22" are adopted. We 
also adopt a terminal velocity projected on the line of 
sight, V cxp =120 km s _1 (e.g, Likkel et al. 1992; Sahai 
et al. 1999; Claussen et al. 2009). For an O-rich circum- 
stellar envelope harboring H 2 and OH maser emission, 
/co = 5 x 10~ 4 is adopted. We assume s(3) =0.43 on the 
basis of previous observations with the CSO 10 m tele- 
scope (Knapp et al. 1998; Groenewegen et al. 1999). 
The unknown factor W (0.1^ W <5) introduces an un- 
certainty of the mass loss rate within a factor of 0.4-2.1. 

Thus we derived the mass loss rate of 116342, M gas w 
2.9 x 1O~ 5 M yr" 1 . If we adopt a HWHM of the 
CO spectrum to be the terminal expansion velocity, 



Kxp =79 km s" 
- 5 M Q yr" 1 , 



then we obtain a value M gas w 1.3 x 
10 -0 Mq yr _i , comparable to that previously estimated 
with CO J =2-1 emission (M gas « 1.7 x 1O _5 M yr _1 ob- 
tained after correction 1 , He et al. 2008). Note that these 
calculations assume spherically symmetric mass loss. As 
mentioned later, because the CO emission may be asso- 
ciated with a collimated, fast outflow, the mass loss rate 
is expected to be overestimated. Note also that the de- 
rived mass loss rate is much lower than roughly estimated 
with mid-IR emission (M gas w 10~ 3 Mq yr -1 , Dijkstra 
et al. 2003). Because the latter adopted a much lower 
expansion velocity (15 km s" 1 ), leading to a higher mass 
loss rate with a higher adopted expansion velocity, there 
exists a large discrepancy between the rate values. Zijlstra 



The previous calculation adopted a beam size of the IRAM 
30 m telescope to be 21", but which should be corrected to be 
12". When the derived mass loss rate was scaled to that of the 
Arizona Radio Observatory 10 m telescope (32"), we obtained 
the new value for the CO J = 2 — 1 emission. 
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et al. (2001) argues the possibility that the value derived 
in the latter is attributed to the extreme extinction com- 
ing from an edge-on thick disk, which may be invisible 
in the CO emission as explained later. Using any calcu- 
lated mass loss rate, it is concluded that 116342 should 
be at the phase of stellar evolution, where its mass loss 
is highest. The success of the CO emission detection to- 
wards 116342 may be attributed to its close proximity to 
the Earth and the large offset from the Galactic plane, 
which enables us to more easily obtain a long integration 
time more easily. Even with a similar mass loss rate and 
a similar source distance, we need a much longer obser- 
vation time towards W 43 A including adjacent off-target 
positions in the Galactic plane. Alternatively, interfero- 
metric observations enable us to obtain higher sensitiv- 
ity as achieved with the Atacama Large Millimeter and 
Submillimcter Array (ALMA) and to distinguish intrinsic 
CO emission from the background emission. 

4-2. Copious mass loss from an expanding flow in 
IRAS 16342-3814 

In the case of an optically thick spherically-expanding 
circumstellar envelope, the observed CO emission shows 
a parabolic spectral profile, which has been confirmed in 
many AGB stars (c.f, Kemper et al. 2003). When a bipo- 
lar outflow is a dominant component of the CO emission, 
the spectrum exhibits double-sided wings whose full veloc- 
ity width indicates twice the velocity of the bipolar flow. 
Such wing profiles are found in the spectra towards young 
stellar objects (YSOs) and the spectra themselves resem- 
ble Gaussian profiles. More intense CO emission around 
the systemic velocity suggests that the bipolar outflow is 
accelerated in the vicinity of a YSO. Even if the YSO is 
deeply embedded in a dense molecular cloud core, the core 
itself may not be detected in the higher J transition of CO 
emission when it is too cold to excite the emission and/or 
its filling factors are too small in the single-dish beam. 
The CO spectrum of II6342 resembles such a Gaussian 
profile, however it is difficult to expect a physical and dy- 
namical condition similar to those of star forming regions. 

In order to better understand the morpho-kincmatical 
structure of the 116342 CO flow, we have constructed mod- 
els using the Shape software package (see, e.g., Steffen & 
Lopez 2006). Shape is a morpho-kinematic modeling tool 
used to create three-dimensional(3D) models of astronom- 
ical nebulae. It has also been used for obtaining models of 
circumstellar CO envelopes (e.g., Nakashima et al. 2009). 
Note that Shape does not calculate the numerical hydrody- 
namic evolution, temperature profile, or radiation transfer 
equations in detail. Instead, a monte-carlo algorithm is 
introduced, in which a three-dimensional space is divided 
into cubes with a specific temperature and contributions 
to emission, absorption, and scattering from the individ- 
ual cubes are calculated. Because Shape has been updated 
since the calculations in Nakashima et al. (2009), we are 
now concerned with the morphology including a density 
profile, kinematics as function of distance from the cen- 
tral star, and opacity, but do not try to reproduce the 
intensity scale. 



Figure 4 shows one of the best-fit models and its simu- 
lated spectrum obtained from an accelerated bipolar flow 
model. The modelled flow is expanding with a velocity 
proportional to the distance from the central star and 
with an exponentially decreasing gas density. The flow 
is slightly squeezed near the central star in order to repro- 
duce the slightly flat peak of the observed spectrum. The 
simulated brightness distribution resembles the shape of 
the Keck optical image (Sahai et al. 2005), except for the 
central brightest part which is obscured in the latter due 
to heavy extinction. Note that it is impossible to obtain a 
unique model to reproduce the observed spectrum without 
any information on observed spatial brightness distribu- 
tion. We confirmed that there exists other possible models 
to reproduce the same Gussian spectral profile, but which 
are based on different flow morphologies and density pro- 
files. 

Nevertheless, through our Shape modelling, we find pos- 
sible major factors to control the spectral profile and can 
qualitatively evaluate them. First, as mentioned above, 
the CO emission should be optically thin. The II6342 CO 
outflow cannot be approximated to a spherically symmet- 
ric flow. Taking into account the flow major axis that may 
be coincident with the major axis of the H2O maser jet 
at an inclination of 44 with respect to the line of sight 
and at a position angle of 66 east from north (Claussen 
et al. 2009), an asymmetric profile is easily formed if the 
CO emission is optically thick. Also taking into account 
the 3-D speed of 180 km s _1 and the acceleration (or 
the velocity gradient) of the jet, an optically thin flow is 
more reasonabl. Secondly, the CO flow should be domi- 
nated by a bipoloar flow rather than an expanding torus 
or an equatorial flow. Nakashima et al. (2009) repro- 
duced the spectrum and spatial brightness distribution of 
the post-AGB star IRAS 07134+1005 in term of a torus 
expanding with a velocity of 8 km s _1 . When adopt- 
ing a similar morpho-kinematical structure, we find that 
the simulated spectrum has too high brightness at the 
high velocity wings. Verhoelst et al. (2009) proposed 
the existence of a dusty equatorial flow that has a bicon- 
ical cavity along the H2O maser jet axis and a full flow 
opening angle of 145°. Assuming that the CO emission 
originates from the same equatorial flow, we fail to re- 
produce the Gaussian wings even by increasing the flow 
terminal velocity up to 150 km s . Instead, by assuming 
the CO emission originating from the cavity of the equa- 
torial flow, we succeed in reproducing the Gaussian wings. 
This model is also consistent with the velocity distribution 
found in 1612 MHz OH maser emission, in which the high 
velocity maser components are significantly spatially de- 
viated from the collimated jet found in H2O maser emis- 
sion (Sahai et al. 1999; Zijlstra et al. 2001). Thirdly, 
a dense torus or an equatorial flow proposed by Dijkstra 
et al. (2003) and Verhoelst et al. (2009), respectively, 
may contribute to the emission or self-absorption of the 
J = 3 — 2 CO line from 116342, if it exists. In the equatorial 
flow model mentioned above, instead of reproducing the 
Gaussian wings as observed, a bump at the spectrum peak 
appears. It may be decreased by self-absorption, but this 
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causes an unwanted asymmetry in the profile. A decrease 
in gas density near the central star also reduces the bump, 
but this is inconsistent with the existence of a thick torus 
that is clearly seen in optical and mid-infrared images as 
a dark lane (Sahai et al. 1999; Dijkstra et al. 2003; Sahai 
et al. 2005; Verhoelst et al. 2009). Here we suppose that, 
as mentioned above, the CO outflow is slightly squeezed 
to reduce the contribution of the central part of the mod- 
elled CO spectrum (figure 4a) . The negligible contribution 
from the dense torus to the observed CO spectrum may be 
attributed, as mentioned above, to the filling factor of the 
torus that may be much smaller than that of the CO out- 
flow and/or to the too cold (^20 K). The large difference 
in the filling factors is expected if the total length of the 
CO flow is larger than that of the H2O maser jet and the 
optical lobes (~3", Sahai et al. 2005; Claussen et al. 2009) 
and/or if the radius of the equatorial flow is limited within 
1" (Verhoelst et al. 2009). For the temperature structure, 
He et al. (2008) demonstrated the existence of warmed 
walls of the bipolar lobes for generating a wing component 
of CO. Verhoelst et al. (2009) also confirmed the extreme 
extinction should be attributed to cooler dust around the 
equatorial flow. Note that, the mass loss rate caluca- 
tion in equation (1) assumes a constant expansion 
velocity of the flow while the model adopts a linear 
acceleration of the flow. Even for such a linearly 
accelerating flow, a mass loss rate can be defined 
if the density profile has a form of pcxr 3 . In the 
model, the outflow density decreases exponentially 
with distance, leading to a good approximation to 
the r 3 profile. 

Thus, we have obtained several insights into the high 
velocity CO outflow in 116342 through the SHAPE mod- 
elling. When we obtain radio interfcrometric images as 
obtained with the SMA and ALMA, we will be able to 
obtain a unique morpho-kinematical model with more pre- 
cise information on the gas density structure to completely 
explain the observed outflow. This should also be the 
first step to understanding the process of the (almost) si- 
multaneous developments of spherically-expanding, slow 
tori and highly-collimatcd, fast jets as well as their de- 
cays in the later evolutionary phase. It has been demon- 
strated that a developed torus as seen in the 116342 dark 
lane should be created with a highly collimated jet ei- 
ther almost simultaneously or before the jet launch with 
a very short time lag (^300 yr, Huggins 2007). In fact, for 
116342, dynamical time scales of the torus and the jet are 
estimated to be 320 yr (1000 AU/15 km s" 1 ) and 110 yr 
(3000 AU/120 km s^ 1 ), respectively. The evolution of 
the expanding torus should be explored on the basis of a 
larger number of sampled sources in different evolutionary 
phases. 116342 should be an important sample for statis- 
tical analyses as demonstrated by Huggins (2007). Some 
discrete mass ejection events such as companion interac- 
tions arc also expected in order to explain the simultane- 
ous evolution of jet and torus. It is also worthwile to note 
that the CO emission in 116342 originates from a different 
component than that of IRAS 07134+1005, which may 
be in a later evolutionary status than 116342 in the post- 
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AGB phase (Nakashima et al. 2009). Not only water 
fountain sources but also H2O maser sources associated 
with planetary nebulae should be the targets for the ra- 
dio interfcrometric observations in order to discriminate 
possible scenarios of the evolution and de-evolution of the 
stellar jets and the circumstellar tori. 

5. Conclusions 

Using ASTE, we have surveyed CO J = 3 — 2 emission 
from nine water fountain sources and detected the CO 
emission towards 116342 and IRAS 18286-0959. 116342 
may have a mass loss rate of M gas « 5.6 x 10~ 5 M yr -1 , 
with most emission attributed to a bipolar outflow rather 
than a high velocity expanding torus or a biconical flow 
with a wide opening angle. On the other hand, to further 
clarify the intrinsic CO detection from IRAS 18286-0959 
with a very narrow velocity width (<~ 3 km s" 1 ), we 
need interfcrometric observations of this object. Because 
the water fountain sources are located at large distances, 
higher sensitivity is essential in future observations. 
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Table 1. Parameters of the water fountain sources and the ASTE observations. 







l L 






D 4 


tjet 




Dur' 


On- 8 


rms 


IRAS name 


Other name 


(") 


(km s 


- 1 ) 


(kpc) 


(year) 


Ref. 6 


(hr) 


point 


(mK) 


15445-5449 


OH 326.5-0.4 




- -100 a 


111 


-7 


? 


2 


0.90 


1 


10 


15544-5332 


OH 325.8-0.3 




--74 9 


40 


3.5 or 10.4 11 


? 


2 


0.90 


5 


62 


16342-3814 


OH 344.1+5.8 


2.4 


42+2 


240 


2.0 


-100 


1 


4.15 


1 


4 


16552-3050 


GLMP 498 


0.11 


-12 10 


170 


19.6 11 


-60 


8, 9 


0.75 


1 


8 


18286-0959 




0.24 


41+1 


200 


3.1 


-15 


3, 4 


2.10 


5, 9 


35 


18450-0148 


W43A, OH31.0+0.0 


0.8 


34+1 


190 


2.6 


50 


6, 7 


1.10 


5 


20 


18460-0151 


OH 31.0-0.2 


0.11 


125+1 


290 


2.0 


-5 


3, 4 


1.00 


9, 1 


34 


18596+0315 


OH 37.1-0.8 




90+1 


59 


4.6 or 8.8 11 


? 


2 


0.90 


1 


19 


19134+2131 




0.15 


-65+2 


100 


8.0 


40 


5 


1.50 


1 


4 



1 Total angular length of the jet system. 
2 Systemic velocity of the jet system. 

3 Full range of the linc-of-sight velocities of water maser emission. 

4 Distance to the source. 

5 Dynamical age of the jet (« l/AVi os ). 

6 References of the jet parameters. 1: Claussen et al. (2009); 2: Deacon et al. (2007); 3: Deguchi ct al. (2007); 4: 
Imai (2007); 5: Imai et al. (2007); 6: Imai et al. (2005); 7: Imai ct al. (2002); 8: Suarez et al. 2007; 9: Suarez 
et al. (2008) 

7 Duration of the total observation time with ASTE. 

8 Number of points observed on and around the target, except an off-point. 

9 Both of the OH and H 2 maser spectra have either irregular or single-peaked profiles, causing difficulty in estimating 
a precise systemic velocity. 

10 No OH maser is detected and the H2O maser spectrum has an irregular profile, causing difficulty in estimating a 
precise systemic velocity. 

11 Kincmatic distance estimated on the basis of the Galactic kinematical parameters provided by Reid et al. (2009) 
and Sofue et al. (2009) is given. 
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Fig. 1. CO J = 3-2 omission spectrum of IRAS 16342-3814. 
The horizontal grey dashed line shows the zero-tcmpcraturc 
baseline. The dotted line shows the Gaussian fit with a peak 



of T, 



MB 



: 32.9 ± 1.0 mK at LSR velocity of 50±2 km s" 1 ^ 



vertical dashed line) and a HWHM width of 79±3 km s 
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Fig. 3. Spectra of CO J = 1 — emission towards 
IRAS 18286—0959 and the adjacent points observed with 
the Nobeyama 45 m telescope. (a) Raw spectra at 
IRAS 18286—0959 and four adjacent points, whose relative 
positions are described at the right sides of the spectra. The 
antenna temperature values of the spectra are shifted in a step 
of 3 K for clarity. There exists a narrow peak at Vlsr— 65 
km s~ x , which is brightest at (0", 0") and expected to be 
the intrinsic emission from IRAS 18286—0959. Two wide 
bumps at Vlsr— 50 km s" 1 and V^gR— 90 km s _1 may 
be the Galactic background CO emission, (b) Spectrum to- 
wards IRAS 18286—0959 obtained after an average spectrum 
synthesized from the adjacent points is subtracted. 



No. 1 CO J = 3 — 2 emission from the water fountain sources 
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Fig. 2. CO J = 3 — 2 emission in IRAS 18286—0959. (a) Temperature variations of CO emission peaks along the position offset 
from IRAS 18286—0959 in the east-west (X) direction. The dashed line shows the Gaussian fit for the 65 km s peak with a 
peak of T MB = 1.58 ±0.03 K at a peak offset X =4 f'3±0 '.'2 and a HWHM width of 27 !'3±0 '.'6. (b) Same as (a) but in the 
north-south (Y) direction. The dashed line shows the Gaussian fit for the 55 km s — 1 peak with a peak of Tjyjg = 0.79.8 ±0.33 K 
at a peak offset of Y =5 ."1±3 ."8 and a HWHM width of 10 ."6±7 ."3. (c) Spectrum at the position of IRAS 18286-0959, which is 
subtracted by the mean spectrum synthesized from spectra at four positions around IRAS 18286—0959. (d) Same as (c) but enlarged 
in velocity around the emission peak. The dashed line shows the Gaussian fitting slope with a baseline offset of 0.04±0.01 K, a peak 
of T MB = 1.07±0.06 K at LSR velocity of 65.3±0.1 km s~\ and a FWHM width of 3.0±0.2 km s _1 . 
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(C) 



Fig. 4. A SHAPE model that reproduces the observed CO 
spectrum, (a) Mesh image of the modelled morphology, (b) 
Brightness images of the model, (c) Spectrum obtained from 
the model (a black line), whose velocity width and peak an- 
tenna temperature are scaled to those of the observed spec- 
trum (a grey line). 



